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Abstract 

We study a smooth crossover transition in the gas of quark- gluon bags. The equation of state at 
high temperature is that of the quark-gluon plasma. However, the system consists of the bags with 
finite volumes which are defined by the model parameters of the mass- volume bag spectrum. Possible 
structures in this cluster quark-gluon plasma are classified. 
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I. INTRODUCTION 



A connection of the bag model [1| with statistical description of strongly interacting matter 
at high energy density has a long history [2]. The possibility of phase transitions in the gas 
of quark-gluon bags was demonstrated for the first time^in^Ref.^ 3). Further studies allowed 
to obtain the 1st, 2nd, and hig 



transitions was also pointed out 



rer order transitions 

5|. Recently it was suggested |9| to model a smooth crossover 
transition by the gas of quark-gluon bags. Inspired by this suggestion we study in more details 
the high temperature behavior of the system of quark-gluon bags in case of the crossover. Note 
that the present and future experimental facilities such as RHIC and LHC produce strongly 
interacting matter in the crossover region of the QCD phase diagram |l0). Moreover, the 
ocation of the (tri)critical point that ends the line of the 1st order phase transition is unknown 



. It can happen that nucleus-nucleus collisions at SPS and FAIR also partially enter the 



crossover region. 

Section [TT] gives a short overview of the gas of bags with excluded volume. In Section [TTT1 the 
phase diagram of the model parameters in the mass-volume bag spectrum is considered. The 
average volume of the bags is calculated in the high temperature limit. In Section [IV] a short 
summary of the results is presented. 



II. GAS OF QUARK-GLUON BAGS 

The partition function for gas of quark-gluon bags is the following j^]: 

oo N / N \ N / N \ 

Z(V,T) = m II / dm * dv * -°( m *>0 <^ T ' m *) [ V - 9 [ V - 5>i 

N=0 i=l "* \ j=l J V j=l J 



(1) 



where V and T are the system volume and temperature respectively, N is number of bags, m,j, 
Vi are mass and proper volume of i-th bag, p(rrii, Vi) is the mass-volume spectrum of bags that 
will be specified later. The 0(T, m,) in Eq. (TjQ) equals: 



4>(T,mi) = 



2tt 2 



k 2 dk exp 



(k 2 + m 2 ) 1 ' 2 



T 



m 2 T frrii 



(2) 



and it has the physical meaning of the particle number density in Boltzmann ideal gas, K2 in 
Eq. (T5]) is the modified Bessel function. The baryonic number and other conserved charges are 
assumed to be equal to zero. The equation of state can be most easily studied with the help of 
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the Laplace transform: 

/>oo 

Z(T,s) = / dVexp(-sV) Z(V,T) = [s - f(T, s)}^ 1 , (3) 
Jo 

where 

POO 

f(T,s) = / dmdv exp(— vs) p(m, v) 0(T, m) . (4) 



Jo 

In the thermodynamic limit, V — > oo, the partition function behaves as Z(V,T) = 
exp [pV/T], where p(T) is the system pressure. An exponential increase of Z(V,T) on V gen- 
erates the singularity s* of the function Z(T, s) in variable s. Consequently, one can calculate 
the pressure knowing only the position of the farthest-right singularity s*: 

P (T) = T lim hlpll = T s*(T) . (5) 

There is a pole singularity s* = sh of Z(T, s) calculated from the transcendental equation: 

s H {T) = f(T,s H {T)) . (6) 

Another singular point of Z(T, s) denoted as sq(T) emerges due to a singularity of the function 
f(T, s) itself. The system pressure takes then the form 3j: 

p(T) = Ts*{T) = T-max{s H (T),s Q (T)} , (7) 

i.e. the farthest-right singularity s*(T) of Z(T,s) © can be either the pole singularity sh{T) 
([6]) or the sq(T) singularity of the function f(T,s) (H]) itself. The mathematical mechanism 
for possible phase transition (PT) in the gas of quark-gluon bags is the "collision" of the two 
singularities, i.e. sh(T) = sq(T) at the PT temperature T = Tq. 

The crucial ingredient of the model is the form of the mass- volume spectrum p(m, y). In the 
case of a bag filled with the non-interacting massless quarks and gluons 1 one finds 0, Q, S] : 



p(m,v) ~ C v 1 {m — Bv) & exp 



- <t* /4 v 1 / 4 (m - Bvf' 4 
3 



where C, 7, 5 and B, the so-called bag constants, are the model parameters, and oq = 957r 2 /60 
is the Stefan-Boltzmann constant counting gluons (spin, color) and (anti-)quarks (spin, color 



1 This picture is reasonable in the region where both the mass of the bag, m, and the volume of the bag, v, are 
large. A general form of the mass- volume spectrum function should contain also the low-lying hadron states. 
However, the type of the PT or the high temperature behavior when the PT is absent are not sensitive to 
presence of these hadron-like states. 
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and u, d, s-flavor) degrees of freedom inside the bag. This is the asymptotic expression assumed 
to be valid for a sufficiently large volume and mass of a bag, v > Vn and m > Bv + Mq. The 
validity limits can be estimated to be Vo ~ 1 fm 3 and M ~ 2 GeV [5|. 

The integral over mass in Eq. (jlj) can be calculated by the steepest descent estimate and 
one finds: 

' rr\ 3/2 

ml \ 



roc roc 

/(T, s) ~ C / dv v 7 exp [ — v (s — sq)] / dm(m — Bv) 
Jv ' J Mo 



In ) 



cxp 



(m — m) 2 
8ctqvT 5 



poo 

~ u{T) \ dvv 2+1+s exp[-v(s -s Q (T))] , 

Jvo 

where m = v(a Q T 4 + 5), u{T) = Cn- l a+ l/2 T 4+4<5 ((TqT 4 + B) 3/2 and 



(9) 



sq(T) = ^ (T Q T 3 - — . 



(10) 



III. THE 7-5 PHASE DIAGRAM 



The model parameters 7 and 5 define the presence, location and order of the PT in the gas 
of quark-gluon bags. This is illustrated by the 7 — 5 phase diagram in Fig. [T] (left). 




1a v —V 



-7/4 1b _ 

v = const 



1c p— 00 



3c I 




PT / 




II/ 




I ^ 






^ // 





FIG. 1: Left: The 7 — (5 phase diagram for the gas of quark-gluon bags. Right: The schematic view of 
temperature dependence of the average volume of the bag in different regions of 7 — 5 phase diagram. 
See text for details. 



Most interest in the previous studies 



PTs are discussed in details in Ref. 



anaaaii 



was devoted to PTs. Different order 



5fl. We concentrate here on the situations when PTs are 
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absent. Eq. (ED for Sh(T) can be written as follows: 

POO 

SH = U {T) / dvv a ~ x exp(-v As) oc T 10+45 (Asp T(a,V As) , (11) 

where a = 7 + 5 + 3, As(T) = sh — sq, and T(a, 6) is the incomplete Gamma-function. The 
function f(T,s) has the singular point s = sq, and J(T,sq) can be either finite or infinite 
depending on the value of 7 + 5: 

1) 7 + 5 > -3, 2) 7 + 5 = -3, 3) 7 + 5 < -3. (12) 

The phase transition never exists for 7 + 5 > —3 because of f(T, sq) = 00, and in this case 
the solution sh(T) of Eq. ffTTl) always corresponds to s# > Sq, see Fig. [2] left. For 7 + 5 < — 3 
the existence of a PT depends on the value of parameter 5. There are three distinct cases: 

a) 5 > -7/4 , b) 5 = -7/4 , c) 8 < -7/4 . (13) 





s s 

Q H 




s s s = s 

Q H Q H 



FIG. 2: The solid lines present the dependence of / on s at fixed temperatures. The pole singularity 
sh and singularity sq are denoted by crosses and circles, respectively. Left: The cases 1 and 2 in 
Eq. (|12|) correspond to f(T, sq) = 00. Middle: The cases 3a and 3b in Eqs. (|12|13|) are shown. They 
correspond to f(T, sq) > sq at all T. Right: The case 3c with T\ < T2 = Tq < T3. It leads to the 
"collision" of two singularities sh = sq at the PT temperature Tq- 



The PTs take place in the system of quark-gluon bags in the only case 3c, i.e. if 7 + 5 < —3 
and 5 < —7/4. This is the lower left corner of the 7 — S phase diagram marked as 3c and shown 
by grey color in Fig. [1] left. The region with 1st and 2nd or higher order PTs are marked as I 
and II correspondingly. 

The average volume of the quark-gluon bag can be calculated as: 

f dvdmv p(m,v) (f)(T,m) exp(-s*v) ^ 1 T(a + 1, V As(T)) 



viT) 



J dvdm p(m,v) 0(T, m) exp(— s*v) 



As(T) T(a,V As(T)) 



(14) 
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It can be proven that sh ~ T 3 for T — > oo. Using the asymptotic expansion for incomplete 
T-function (see Appendix) one obtains from Eqs. ( TTTT) . (HM the behavior of average volume 
v{T) of the quark-gluon bag at large temperature T: 

la, 2a, 3a : As(T) ~ InT -> oo , v(T) -> V , (15) 

16, 26, 36 : As(T) S const > v{T) -> const , (16) 

lc : As(T) ~ T( ?+4<5 )/ a -> , U(T) ~ T^ 7 ^^ -> oo , (17) 

2c : As(T) ~ exp (- T _7_4<5 ) -> , v(T) ~ exp (7 1 - 7 - 44 ) -> oo . (18) 

The results ( jToTfTBl) for u(T) are schematically presented in Fig. [1] right. Note that if the 
phase transition takes place at T = Tc (case 3c) the v(T) becomes infinite at T > Tq in the 
thermodynamic limit. This is also shown in Fig. [T] right. 



IV. SUMMARY 



In this paper we have studied the gas of quark-gluon bags at high temperature T. The 
mass- volume spectrum function p(m,v) for the quark-gluon bags is taken in the form (jHJ). The 
behavior of the system depends crucially on the values of the 7 and 6 parameters in Eq. (jHJ). 
The special regions of the 7 — 5 phase diagram defined by Eqs. (112j|13p lead to different behavior 
of the gas of the quark-gluon bags. The pressure p and energy density e for different values of 
7 and S have the same asymptotic behavior at high temperature, p = <tqT 4 /3 and e = <JqT a . 
This corresponds to the equation of state of non-interacting quarks and gluons inside the bags, 
i.e. ideal quark-gluon plasma (QGP). However, the average volume of the bag v(T) (1151I18P 
and average mass m(T) = w(ctqT 4 + B) have rather different behavior in different regions of 
the 7 — 5 phase diagram. If the system of quark-gluon bags has no phase transition the v(T) 
remains finite at high temperature. Such a cluster QGP can be rather different from the ideal 
QGP despite of the similar to that equation of state. The kinetic properties of the cluster QGP, 
e.g. the shear and bulk viscosity, may deviate strongly from those in the quark-gluon gas. 
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APPENDIX A 



The incomplete gamma function F(a,x) (see, e.g., Ref. |12|) has the asymptotic expansions 
at x — > oo, 



T(a,x -> oo) = z^e - * 1 + 



.r 



(Al) 



and x — > 0, 



r(a,x^0) = T(a) (1 + O(x)) 



a^-n , (A2) 



( — l) n x~ n 

= - — f-(ifi(n + i) - lnx) H (1 + 0(x)) , a = -n , 

n\ n 

where ip(z) is the logarithmic derivative of the gamma function, ip(z) = T'(z)/T(z). Eq. ( 1A2 
gives for a = 0, —1, —2: 



T(a, x — > 0) 



-mx + 0(1) , 

l + 0(x) 
\nx H — 

1 + Q(x) 
2x 2 ' 



a = 



a = -1 



a = -2 . 



(A3) 
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